Field-induced chiral phase in isotropic frustrated spin chains 
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It is shown that an external magnetic field applied to a spin-5 isotropic zigzag chain induces 
a phase with spontaneously broken parity, characterized by long range ordering of vector chiral- 
ity. This is in contrast to the two-component Luttinger liquid scenario proposed in the literature. 
Relevance to real materials is discussed. 
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Introduction. - In recent years, phases with broken vec- 
tor chirality in quantum spin chains have attracted a con- 
siderable interest 0,0) El El HI- They are characterized by 
a nonzero expectation value of the vector product of two 
adjacent spins 4 = (S n x SVi+i), so that in a chirally or- 
dered phase spins tend to "rotate" in a certain preferred 
plane predominantly clockwise or counterclockwise. This 
type of order breaks only a discrete symmetry between 
left and right and thus is allowed in one dimension (ID), 
in contrast to the long-range helical spin order [(| 0| ■ Re- 
cently, chirally ordered phases were numerically found in 
frustrated chains with easy-plane anisotropy 0, 0] ■ The 
chiral ordering transition has possibly been observed ex- 
perimentally in the quasi-lD anisotropic organic magnet 
Gd(hfac) 3 NITiPr 0]. 

In this Letter we show that a chiral phase emerges 
in isotropic frustrated spin chains as well, if they are 
subject to a strong external magnetic field. We focus on 
the model of a zigzag chain defined by the Hamiltonian: 

n = j 1 J2s j - s j+1 + j 2 J2 Sj ■ s j+2 -hj^s* (i) 

3 3 3 

where Sj are spin-5 operators at the j-th site, and 
J\,2 > 0. It is easy to analyse the classical counter- 
part of the above model, where spins are represented by 
vectors, (S^,S*) h-> (S , sin0„e ±46 '", S'cos^n). The ap- 
plied field selects a preferred plane, reducing the sym- 
metry to U(l). Depending on the frustration strength 
a — J2/J1, the in-plane ground state configuration is 
given by 6 n = (w ± X)n, with A = for a < 1/4 
and A = arccos(l/4a) for a > 1/4, respectively. The 
spins are canted towards the field, cos0„ = h/h Sl where 
h s = AS{ J\ cos 2 (A/2) + J2 sin 2 A} is the saturation field. 
The classical ground state is a canted antiferromagnet 
for a below the Lifshits point |, while for a > j one 
has two degenerate helical ground states, as reflected by 
the ± signs above, which correspond to the left and right 
chirality k = ±5 2 sinA. Thus, for a > j in presence of a 
field the initial SU(2) symmetry is reduced to U(l) x Z2. 

In the quantum case the U(l) symmetry cannot be bro- 
ken, but it is allowed to break the discrete Z2 chiral sym- 
metry. Such a scenario is realized in anisotropic chains, 



where two different chiral phases were found 0,0,0, 0,0]. 
A natural question arises: Can an external magnetic field 
act similarly to the xy anisotropy j^] , favoring the chiral 
order in isotropic spin chains? Recent numerical studies 
|lfl llll | propose the scenario of a two-component Lut- 
tinger liquid without any breaking of the Z2 symmetry, 
casting doubts on the above idea. The aim of the present 
Letter is to show that the correct high-field physics of 
isotropic frustrated chains is indeed determined by the 
spontaneous breaking of the chiral symmetry. 

S — ^ chain.- We start with the extreme quantum 
spin-i case which admits a field-theoretical description 
based on the bosonization approach. Consider the limit 
of strong frustration a 3> 1 and strong magnetic fields 
h ~ J2. The system may be viewed as two chains 
weakly coupled by the zigzag interaction J±. A single 
spin-i chain in uniform magnetic field is known to be 
critical, its low-energy physics being effectively described 
by the standard Gaussian theory |12( known also as the 
Tomonaga-Luttinger liquid: 

U= \\ ^{^(^) 2 + A W) 2 }- (2) 

Here <f) is a compactified scalar bosonic field and is 
its dual, dt<p = vd x 9, with the commutation relations 
[4>{x),0(y)] = iO(y — x), where <d(x) is the Heaviside 
function and the regularization [(f>(x), 0(x)] = i/2 is as- 
sumed. Integrability of the S — | chain model makes 
possible to relate explicitly the coupling constants of the 
theory, the spin wave velocity v and the Luttinger liquid 
(LL) parameter K, to the microscopic parameters J2, 
h. The exact functional dependences v(h) and K(h) are 
known (see |l.^ and references therein) from the numer- 
ical solution of the Bethe ansatz integral equations 
Particularly, K increases with the magnetic field from 
K(h = 0) = 5 to K = 1 for h approaching the saturation 
value 2J 2 . 

In the infrared limit the following representation of the 
lattice spin operators holds [l^ : 

Si = —=d x (f> + — sin \2kpx + V47T0} + m (3) 

S- = {-\) n e - ie ^ {c + b sin (2k F x + V^(j))}, 
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Here m(h) is the ground state magnetization per spin 
which determines the Fermi wave vector fep = (5 — mk 
and is known exactly from the Bethe ansatz results |14j . 
Nonuniversal constants a, b, and c for general h have been 
extracted numerically from the density matrix renormal- 
ization group (DMRG) calculations (lif. 

We treat the J± interchain coupling term perturba- 
tively, representing two decoupled chains in terms of 
Gaussian models of the form (J2J. It is convenient to 
pass to the symmetric and antisymmetric combinations 
of the fields describing the individual chains, <f)± = 
(0i ± 4> 2 )/V2K and 6± = (6>i ± B 2 )^/kJ2. The effec- 
tive Hamiltonian describing low-energy properties of the 
model takes the following form: 



Tin + TLi, 



Ht = -[(9A) 2 + (9^±) 2 ], 



7~Li n t = <7i cos kF cos {JtF + y/8TvK-(f>-^ 

- g 2 d x 9 + sm(y/2Tr/K-e-). (4) 

Only the relevant terms are shown here 0] , including the 
"twist operator" with nonzero conformal spin The 
Fermi velocity v cx J 2 , while the couplings 30,1,2 oc J\ <C 
v. The renormalized LL parameter is given by 



K- = K 



(h){l + JiK(h)/(m>(h))}. (5) 



Note that in the first order of J\j J 2 the correction to K_ 
for the zigzag type of interchain coupling is twice larger 
compared to that for the ladder type of coupling. 

The inter-sector part of contains a term which can 
be identified as an infrared limit of the product of in- 
chain and interchain chiralities: one can show that 



2i,2i+2) K 2i,2i+1; 



(6) 



where k*j = (Si x Sj) z . All the other terms, omitted in 
Q, are made either irrelevant or incomensurate by the 
external magnetic field. The Hamiltonian (0J gives the 
minimal effective field theory describing the low-energy 
dynamics of a strongly frustrated (a » 1) spin-i zigzag 
chain for a nonzero magnetization m. For small m the 
LL parameter K_ ~ |, and the inter-sector g 2 term has 
a higher scaling dimension than the strongly relevant g\ 
term in the antisymmetric sector. In this case the sys- 
tem is in a phase with relevant coupling in antisymmetric 
sector, as discussed for the first time in Ref. yj| (later 
dubbed EO phase 01)- I n contrast to that, at h — all 
terms generated by the zigzag coupling are only marginal. 

When h increases, the chirality product operator © 
can become more relevant than the g\ term controling 
the field 0_; the latter term becomes less relevant with 
the increase of h as well as with the increase of the zigzag 
antiferromagnetic coupling J\. To study this situation, 



one can apply a mean field decoupling procedure to the 
inter-sector term in the spirit of Ref. 1]. At the mean 
field level, the interaction 7i; nt takes the form 



T~(.mf — gi cos kp cos(kp + y/87rKZ(j)-) 
-g 2 d x e + (sm J^9J) - g 2 (d x 9+) sin x f]ft 



(7) 



Remarkably, the mean field Hamiltonian reveals a com- 
petition between the basic and dual field terms of the 
form sin(70_) and sin((5#_) with jS — A-k, exactly the 
value where the Ising quantum phase transition takes 
place 0, To find the critical magnetic field h cr 

which corresponds to this transition, we equate the RG 
masses produced by the operators sin(70_) and sin(<5#_), 
substituting the averages in JJJ by their mean-field val- 
ues found from self-consistency conditions. Doing so, one 
obtains an estimate for the renormalized LL parameter 
A"_ at the transition, which happens to be related to the 
glorious "golden mean" q: 

2K-(h cr )=q=(V5 + l)/2, 

This leads to the following equation for h cr : 



K(h cr ) = |{l 



J%K(h cr ) 

TTv(h cr ) 



(8) 



The fact that K(h) is a monotonically increasing func- 
tion [lil [2fl | implies that the critical field decreases with 
increasing the antiferromagnetic zigzag coupling J\ : 



(dher/dJ!) < for Ji > 0. 



(9) 



Numerically solving Eq. JSJ, one obtains that the max- 
imal value of h cri achieved at J\ — » 0, is approximately 
h cr ~ I.7J2, and the spin wave velocity in this limit is 
still of the order of the bandwidth, v(h cr ) ~ 0.6 J 2l which 
justifies the applicability of bosonization formalism close 
to h cr . Within this approach, there is no indication that 
the chiral phase would be destabilized by a further in- 
crease of the magnetic field, so one may conclude that it 
extends from h cr up to the saturation field h s . 

Recently, influence of strong magnetic fields on a spin- 
i zigzag chain was studied numerically by means of the 
DMRG technique HUGH. The authors of Refs. ESUJ 
explain the presence of cusps in the magnetization curve 
m(h), observed for large a — J2/J1, in terms of the emer- 
gent two- component Luttinger liquid phase. Our find- 
ings suggest an alternative scenario, according to which 
the phase above the upper cusp singularity is still de- 
scribed by a one-component Luttinger liquid, albeit with 
a spontaneously broken left-right symmetry. It can be 
shown that the cusp itself originates from the Ising 
transition at the boundary of the chiral phase, and that 
the phase below the lower cusp singularity is also chiral. 
Thus, comparing our results with the DMRG phase dia- 
gram |lfj , we can conclude that the TL2 phase in Fig. 1 of 
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Ref. [lOj should be identified as a chirally ordered phase. 
For high fields, in the limit J\ <C J2 the stability region 
of this phase expands with increasing J 1; in agreement 
with our result ©■ Near the saturation field, this phase 
extends up to the classical Lifshits point J 2 = \J\, which 
again agrees with our large-5 analysis given below. 

5=1 chain.- One can obtain a bosonized descrip- 
tion of a single (unfrustrated) spin-1 chain in magnetic 
field exceeding the Haldane gap A by accessing the pa- 
rameters K and v either directly from numerical DMRG 
studies |22j or, in a more exotic way, from the exact so- 
lution of the integrable 0(3) nonlinear <7-model (NLSM), 
which itself is believed to provide a proper effective field- 
theoretical description |23j. Then, a zigzag 5 = 1 chain 
in the regime of strong frustration J± <C J 2 can be stud- 
ied along the same lines as done above for the spin-i case, 
i.e., treating the zigzag interaction as a perturbation cou- 
pling two LLs. The LL parameter of a S = 1 chain turns 
out to be increasing from the free fermion value K = 1 
at ft. = A with the further increase of the field, so that 
generally for h > A one has K > 1 [22I I23I ] , which resem- 
bles a ID Bose gas [24|. This fact leads to a considerable 
simplification: coupling two LLs with K > 1 by a zigzag 
interchain coupling yields the same effective field theory 
®, but since K > 1 the only relevant term is the product 
of chiralities ©■ Thus, in contrast to the S = 75 case, a 
strongly frustrated spin-1 chain immediately enters a chi- 
ral phase, described by a one component Luttinger liquid 
with spontaneously broken chiral symmetry, as far as the 
external field becomes higher than the gap value. Simi- 
lar NLSM-based analysis readily applies to any integer-,? 
zigzag chain in the limit of a > 1 as soon as magnetic 
field closes the Haldane gap. 

Large-S frustrated chain close to the saturation field. - 
In the vicinity of the saturation field the emergence of chi- 
rality can be analyzed for an arbitrary spin value S. In 
the coherent state path integral representation the effec- 
tive Lagrangian is given by £ = — S^2 n (l — cos(f> n )dtO n — 
(H). Using the ansatz 



(-lrsmo^y 



(10) 



one can pass to the continuum limit, treating ip\,2 as 
smooth fields and keeping only non-oscillating terms. For 
h close to h s the densities of magnons with momenta 
around tt ± A are small, iV'ol 2 *C 1, so one may neglect 
any terms of a higher than quartic order. After rescaling 
the bosonic fields (25) 1 / 2 -0i.2 — ► ipi/21 one arrives at the 
Lagrangian of the form 

c = I dxY, {*«^-^|3*<M 2 + mNv| 2 } 

■* (7=1,2 ""' 

- \ I dxMi^ + i^r+wi^i 2 ^ 2 }, (11) 

which was recently discussed in the context of ID binary 
Bose-condensate mixtures |25j . The Lagrangian param- 



eters are in our case given by 

fj,= h s -h, mT 1 = 8J 2 5sm 2 A, (12) 
u = h s /S, w = 2{u~ 4Ji(l + Jl/Jl) sin 2 A}. 

In the harmonic fluid approach j2(| the field operators 
and densities can be expressed through scalar bosonic 
fields i?, V as |^ CT | 2 = {p a + d x tp a /ir} J2 m e 2im ^ p " x+ip ^ 
and r/v = {p a + d x if a /-K} l / 2 e^" J2 m e 2im ^P° x+ f°\ and 
for /i > the Lagrangian i|ll|) describes two LLs of the 
form J21, with a density-density interaction. In contrast 
to Ref. [25j , in our case the total particle numbers of the 
components ni^2 = f dx\ijji,2\ 2 (which are separately con- 
served) are not fixed, but are chosen by the system so as 
to minimize the energy at /i < 0. It is easy to show that 
for w > the system is unstable against any perturbation 
making pi =/= p 2 : indeed, e.g., for p\ > p 2 the interaction 
term leads to renormalization p CT 1 — s- p CT + {d XT o a )/iT with 
(9 x tpi) > (d x ip2)- As a result, the chiral Z2 symmetry 
breaks spontaneously and one of the bands a gets de- 
pleted. The effective theory is a single Luttinger liquid 
with the parameter K > 1 depending on the dimension- 
less coupling constant 



7 



Po 



2S 



— ( 

sin A V 



4J 2 5(1 - h/h s ) 



1/2 



(13) 



where p = (2pm) 1 / 2 /n is the equilibrium density for 
small p (i.e., in the vicinity of the saturation field). For 
h — ► h s , when 7> 1, the LL parameter tends to 1 and is 
given by K ~ 1 + 4/7, and for 7 <C 1 (which, despite the 
condition po <C L, is formally possible for large S) one 
has K ~ 71-/^7 0]. 

The chirality order parameter is directly related to the 
density difference, k ~ (l^il 2 — |^>2| 2 ) sin A. Neglecting 
the depleted field and using a known expression for the 
density correlator [24|, one obtains the leading asymp- 
totics of the chirality correlation function: 



(k(x)k(0)) ~ ^{ 



h 2K sin^ A 



J2S 



}• (14) 



The longitudinal spin correlator (S z (x)S z (0)) is also re- 
lated to the density and behaves similarly to (|14|) . The 
leading part of the transversal spin correlator can be ex- 
pressed through (ip> (x) tp(0)) and is given by 



(S + {0)S-(x)) ~2Sp ( 



K 
\irp Q x 



2K e i\x 



(15) 



Discussion. - In summary, we have shown that a suffi- 
ciently strong magnetic field applied to a spin-5 isotropic 
J1-J2 zigzag chain induces a phase with spontaneously 
broken Z2 symmetry, which is characterized by the long 
range vector chirality order and emerges immediately be- 
low the saturation field if the frustration strength J 2 / J\ 
exceeds the classical Lifshits point value \. This chiral 
phase is gapless and its low-energy physics is effectively 



4 



described by a one-component Luttinger liquid. Our re- 
sults refute the two-component LL scenario proposed in 
Refs. dSEl, and in fact may necessitate reconsidering 
the phase diagrams of other frustrated spin models, par- 
ticularly of a biquadratic-bilinear spin-1 chain in mag- 
netic field [5^. One may suppose that the SI phase in 
Fig. 2 of Ref. [27| should have a broken Z2 symmetry in 
the region beyond the Lifshits point. To clarify this issue 
we suggest to measure directly the chirality correlator 
(kqK^} in the limit n — > 00 above the cusp singularity. 
For spin-i chain, such a correlator was calculated only 
for very short distances and indicated emergence of 
at least short-range chirality correlations for h directly 
below the saturation field h s . 

The chiral phase should be able to survive finite tem- 
perature effects since it involves breaking of the discrete 
Z2 symmetry. Less trivially, it has also a chance to sur- 
vive the three-dimensional interaction without tranform- 
ing into a usual helical long-range order: as noted by 
Villain 6] , at finite temperatures the chirality correlation 
length is much larger than the spin correlation length, so 
with decreasing temperature the chiral order should set 
in before the helical spin order does. 

Several materials are known which realize zigzag spin- 
i chains (see Table 1 of Ref. [HI)- A promising candi- 
date substance for detecting the field-induced chirality 
would be (N 2 H 5 )CuCl3, since its small exchange con- 
stants Ji ~ 4 K and J2 ~ 16 K make feasible the task of 
attaining magnetic fields comparable to J 2- Experimen- 
tally, the projection of vector chirality k on the applied 
field direction could be detected by comparing inelastic 
scattering intensities for oppositely polarized neutrons, 
as it was done for the triangular lattice antiferromag- 
net CsMnBr3 |30| ; a similar route can be employed with 
polarized light. We hope that our results will stimulate 
further experimental work in this direction. 
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